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Abstract 

In 1959 Buchdahl [13] obtained the inequality 2M/R < 8/9 under 
^SJ ■ the assumptions that the energy density is non-increasing outwards and 

' that the pressure is isotropic. Here M is the ADM mass and R the 

area radius of the boundary of the static body. The assumptions used 
to derive the Buchdahl inequality are very restrictive and e.g. neither 



of them hold in a simple soap bubble. In this work we remove both of 



these assumptions and consider any static solution of the spherically 



5h ■ symmetric Einstein equations for which the energy density p > 0, and 

the radial- and tangential pressures p > and pr, satisfy p + 2px < 
Op, O > 0, and we show that 

2m(r) (l + 20) 2 -l 

where m is the quasi-local mass, so that in particular M = m(R). 
We also show that the inequality is sharp. Note that when = 1 
the original bound by Buchdahl is recovered. The assumptions on the 
matter model are very general and in particular any model with p > 
which satisfies the dominant energy condition satisfies the hypotheses 
with = 3. 

1 Introduction 

The metric of a static spherically symmetric spacetime takes the following 
form in Schwarzschild coordinates 

ds 2 = - e Mr) df 2 + e 2A(r) dr 2 + ^2 + ^2 
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where r > 0, 9 E [0, ir], <p € [0,27r]. Asymptotic flatness is expressed by the 
boundary conditions 

lim A(r) = lim p{r) = 0, 

r~ >oo r— >oo 

and a regular centre requires A(0) = 0. The Einstein equations read 

e- 2X (2r\ r - 1) + 1 = 8irr 2 p, (1.1) 
e" 2A (2r/i r + 1) - 1 = 8irr 2 p, (1.2) 

p rr + (/i r - X r )(p r + -) = 8irp T e 2X . (1.3) 
r 

Here p is the energy density, p the radial pressure and pt is the tangential 
pressure. If the pressure is isotropic, i.e., p = Pt-, a solution will satisfy the 
well-known Tolman-Oppenheimer-Volkov equation for equilibrium 

Pr = ~MP + P)- i lA ) 

In the case of non-isotropic pressure this equation generalizes to 

2 

Pr = -fJrip + p) (p-Pr)- (1-5) 

r 

Note that the radial pressure p is monotone in the isotropic case if p + p > 
since p, r > 0, cf. (|2.2j) . The quasi-local mass m = m(r) is given by 

pr 

m(r) = / 4:7rrj 2 p(rj)dr], (1-6) 
Jo 

and the ADM mass of a steady state for which the energy density has support 
in [0, R] is thus given by M = m(R). 

Schwarzschild asked already in 1916 the question: How large can 2M/R 
possibly be? He gave the answer 2M/R < 8/9 [25] in the special case of 
the Schwarzschild interior solution which has constant energy density and 
isotropic pressure. In 1959 Buchdahl [13] extended his result to isotropic 
solutions for which the energy density is non-increasing outwards and he 
showed that also in this case 



2M/R < 8/9. 

This is called the Buchdahl inequality and is included in most text books 
on general relativity in connection with the discussion of the interior solution 
by Schwarzschild, cf. e.g. [26] and [27]. The quantity 2m /r is fundamental 
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for determining the spacetime geometry of a static spherically symmetric 
spacetime, cf. equations ()2.1|) and (|2,4p . A bound on 2M/R has also an 
immediate observational consequence since it limits the possible red shifts 
of spherically symmetric static objects. 

The assumptions made by Buchdahl are extremely restrictive as pointed 
out by Guven and O Murchadha [T7], e.g. neither of the assumptions hold in 
a simple soap bubble and they do not approximate any known topologically 
stable field configuration. Moreover, astrophysical models of stars are not 
unusually anisotropic. Lemaitre proposed a model of an anisotropic star 
already in 1933 [19], and Binney and Tremaine [9] explicitly allow for an 
anisotropy coefficient (cf. also [H] and the references therein). 

One motivation for this study has its roots in the numerical investigation 
of the spherically symmetric Einstein- Vlasov (ssEV) system J7] which admits 
a very rich class of static solutions. The overwhelming number of these 
have neither an isotropic pressure nor a non-increasing energy density, but 
nevertheless 2M/R is always found to be less than 8/9, cf. [7]. There are 
sometimes arguments which claim that the monotonicity of p is necessary 
for the stability of a steady state, cf. e.g. |27j . but at least for Vlasov matter 
this is not the case by the results presented in [6]. 



In this work the problem of finding a sharp bound on 2m /r is solved in 
full generality in the class of matter models which satisfy 



for any static solution of the spherically symmetric Einstein equations which 
satisfies (II. 7p . The class of matter models defined by (II. 7h is very general. 
Indeed, a realistic matter model should satisfy the dominant energy con- 
dition (DEC) which implies that p > and that the inequality (11. 7p holds 
with = 3. The remaining condition that p is non-negative is a standard as- 
sumption for most matter models in astrophysics. Moreover, Vlasov matter 
satisfies the conditions in (jl.7p with f2 = 1. An interesting feature of Vlasov 
matter, in comparison with a fluid model, is that no equation of state which 
relates the pressure and the energy density has to be specified. For Vlasov 
matter, p,p and pt are all determined by a single density function on phase 
space, cf. [5] and |24| for more information on Vlasov matter and the EV 
system. 



p + 2px < Qp, where f2,p and p are non-negative. 



(1.7) 



We will show that 



sup 

r>0 



2m(r) (1 + 20) 2 - 1 
r - (1 + 20) 2 
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The bound that we obtain for 2m jr is sharp in the sense that an in- 
finitely thin shell of matter, with 2m jr equal to the critical value, will sat- 
isfy a form of the generalized TOV equation which allows p and pr to be 
measures (p = here). This is described in detail in the next section. It 
should here be pointed out that for the ssEV system the results in [3J show 
that there exist regular static solutions with the property that 2M/R takes 
values arbitrary close to 8/9 (Q = 1 for Vlasov matter). These solutions 
do approach the infinitely thin shell mentioned above as 2M/R — > 8/9. In 
section 4 we give an analogy with a classical problem in electrostatics (or 
equivalently in Newtonian gravity) which shares the property that the max- 
imizer is a measure at the boundary. In the work by Buchdahl, the solution 
that maximizes 2M/R is the Schwarzschild interior solution with constant 
energy density. This solution has the property that the pressure becomes 
unbounded as 2M/R — ► 8/9, and therefore the solution does not satisfy the 
DEC and is not a realistic steady state. 

Before finishing this section with a review of previous results, let us point 
out that the original motivation for investigating the Buchdahl inequality in 
full generality comes from its possible role in understanding the formation 
of trapped surfaces. Christodoulou has obtained conditions which guarantee 
the formation of trapped surfaces in the case of a scalar field |14j . and this 
result is crucial for his proof of the weak- and strong cosmic censorship 
conjectures |15| . For more information on this see the introduction in [4j. 

General investigations of the Buchdahl inequality have previously been 
undertaken by Baumgarte and Rendall [8] and Mars, Merce Martin-Prats 
and Senovilla [21]. These studies concern very general matter models and 
they obtain the bound 2m /r < 1. This bound gives little information on the 
spacetime geometry since A —* oo as 2m /r — * 1, and in particular it gives no 
bound on the red shift of a static body. In [3] shells supported in [Rq,R%] 
are considered and it is shown that if the support is narrow then a Buchdahl 
inequality holds (i.e. 2M/R < 1 — e, e > 0). This result is superseded by 
the result presented here but some of the ideas in [3] play an essential role 
in this work. Guven and O Murchadha consider the general case in |17j and 
obtain a bound on 2m/ r in terms of the ratio of the tangential- and the radial 
pressure, which they denote by 7. Their bound degenerates (i.e., 2m/r — > 1) 
as 7 —* 00. (Cf. also [10] for a similar analysis which includes a cosmological 
constant.) It is interesting to note that 7 = 00 for the maximizing solution in 
our work since p = and pr is a Dirac measure at the boundary. Also note 
that 7 — > 00 for the sequence constructed in [JJ, which in the limit gives an 
infinitely thin shell with p = and 2px = p- In this context we mention the 
work [16] where an infinitely thin shell is studied. They obtain the bound 
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2M/R < 24/25. Note that this value agrees with our bound when = 2. 
This is not surprising since their infinitely thin shell satisfies the DEC and 
has p = which in our terminology means that = 2. A similar study is 
carried out by Bondi [12] in the case = 1. 

Furthermore, Bondi [TT] investigates (non-rigorously) isotropic solutions 
which are allowed to have a non-monotonic energy density. He considers 
models for which p > 0, p > p, or p > 3p, and obtains bounds on 2M/R 
strictly less than one in the respectively cases. The isotropic condition is 
however crucial since these bounds are violated for strongly non-isotropic 
solutions as this work shows (cf. also 0], [16] and |12j). 

The paper is organized as follows. In the next section we derive our 
basic inequality which only involves p and we formulate our main results. 
The main ideas of the paper are presented in section 3. In section 4 an 
electrostatic analogy (or equivalently a Newtonian analogy) is discussed and 
the proofs of the theorems are given in section 5. 

2 Set up and main results 

Let us collect a couple of facts concerning the system (|l.ip - (|1.3p . A conse- 
quence of equation (jl.ip is that 



and from (jl.2l) it then follows that 



Adding (jl.ip and (jl.2p and using the boundary conditions at r = oo gives 



In particular if R is the outer radius of support of the matter then 



Pr = {—+ 4:Trrp)e 



(2.2) 




(2.3) 



when r > R. Hence 




(2.4) 
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The generalized Tolman-Oppenheimer-Volkov equation (jl,5p implies that a 
solution satisfies 

(m + 4vrr 3 p)e M+A = / 4tt7 ? V +A (p + p + 2p T )dr]. (2.5) 
Indeed, let S = (m + 4vrr 3 p)e' 1+A . Using ([T3]) and ([23]) we get 

— = 4vrr 2 (p + p + 2p T )e fi+X , 
dr 

and the claim follows since S(0) = 0. 

Let us fix r > 0. Consider (|2.5|) . using (|2.3p we get 

(m + 4vrr 3 p)e- £" WP+J>)e 2A */ 

= e -r^(P+P)^ f ^^ p+p)e2Xd(T (p + p + 2p T )dri, 
Jo 

and we have 

m + 4irr 3 p = / 47TT/ 2 e~ ^ 47ro-(p+p) e da^_^^_^_ 2p T )dr\. 
Jo 

Since p is non-negative we obtain the inequality 

m(r) < / A-KT] 2 e~ h ^ a (P + P) e da (p + p + 2pT)drj. 
Jo 

Using again the non-negativity of p and the inequality (|1,7|) we obtain 

m(r) < (1 + «) / r 47rr/ 2 pe - £ ^^^dq. (2.6) 
•/ o 

Note that only p, and not p and g, appears in this inequality in view of (|2.ip . 
This is our fundamental inequality. 

Let £> be the Borel a— algebra of R + and let M. denote the space of 
non-negative a— finite measures on B such that 2m{r)/r < 1, where m{r) = 
J[o r] dh(rj). Let R > and define the operator Fr : M. — > -R+ by 

_ p dfe(>7) 

F R (/i) = / e H 1 — f^Jdhtr). (2.7) 

y[o,ii] 

With abuse of notation it will be understood that Fr(u), where u is a func- 
tion, is the value obtained by applying F to the measure v where dv = udr. 
Now let p = 47rr 2 /j, and note that the inequality (|2.6p can be written 

m(r) < (l + n)F r (p). (2.8) 
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Furthermore, note that by taking p = and 2pT = p the inequalities 
above become equalities and we can for this special class of solutions define 
a form of the generalized TOV equation which is valid whenever 4irr 2 p = 
heM, 

m(r) = (l + £l)F r (h). (2.9) 

This form of the TOV equation will be used to see that the infinitely thin 
shell which maximizes 2m / r satisfies the TOV equation in the sense of mea- 
sures. 

By a steady state we mean a solution of the Einstein equations (|1 . X [) - 
(jl.3h such that p,p and px are C 1 functions on [0, oo). A steady state of 
course satisfies the generalized Tolman-Oppenheimer-Volkov equation. For 
our purposes it is sufficient that the triplet (p,p,pr) satisfies the integrated 
form (|2.5p of the generalized TOV equation. We say that (p,p,pr) is an 
admissible triplet if: each of these functions is in Ll ([0, oo); Airr 2 ), where 
4-7rr 2 is the weight, equation (12. 5p is satisfied a.e., and there is an Q > such 
that (jl.7p holds a.e. The following theorem is our main result. 

Theorem 1 Consider any admissible triplet (p,p,pr)- Then 

2m(r) (l + 2f}) 2 -l 

sup ( -l<^-L 2.10 

r>o r ~ (1 + 20) 2 v ; 

The arguments leading to Theorem 2 in [3] (and also the arguments in the 
proof of Theorem 1 above) imply that the bound (|2.1U|) is sharp in the sense 
given by the theorem below. Before stating this theorem let us introduce 
the notation for the Dirac measure at r = R. 

Theorem 2 Take R > 0, and let 

i? ((l + 2Q) 2 -l) 
2 (1 + 2Q) 2 ' 

Let 

Mu R 
P ~ i^R 2 "' 

and let p = and 2px = Qp, then 112. 9\) holds with h = AirR 2 p and r = R. 

3 Main ideas 

The details of the proofs make the main ideas become less transparent so 
let us describe them in this section. 
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Given a steady state with support in [0, J?], R > 0, there is a smallest 
r* £ [0, with the property that 

2m(r*) 2m(r) 
= sup ■ 



r>0 r 



We will show that if 



then 



n (1 + 20)5 > 

F - (p)< TTo- (3 - 2) 

In view of (|2,8p we thus obtain a contradiction and no steady state with the 
property (|3.ip can exist. To show that (13. If) implies (|3.2p is of course the 
main difficulty 

We will approximate the given steady state with a sum of step functions. 
The precise way this is done is left to the proof. Let r* be as above and let 

«W = Xft.n]- -j +X[r' v r 2 )— — 7 + "• + X[r'„_ 1 ,r N ]——i > ( 3 - 3 ) 

' 1 ' ' 2 ' i ' N ' N—l 

where {ro, n, rjv} is a sub-division of the interval [Rq, R±], so that ro = Ro 
and tat = r*, and where < rjj. < rk+i, and % is the characteristic function. 
First we take rl = and choose the constants Cj, j = 1,2,... so that u 
approximates p in sup norm. We will then admit the parameters ri to vary. 
Note that 

m u (r k ) _ ]_sr^ 

3=0 

independently of the choices of r'-, where m u (r) := J^udr. 

First we consider the first two terms in (|3.3p and perform the limit 
7q — ► ri and — > r 2 so that the first two step functions become Dirac 
measures at r = r\ and r = r 2 . We then show that the operator F applied 
to the new measure is greater than F(u). More precisely we show that 

F r M<FrM), (3.4) 

where 

u 2 = ciri4-i + c 2 r 2 z5 r2 + ii 2 dr, 

and 

( \ C 3 r 3 , c 4 r 4 , . C N r N 

Mr) = XK ^ ] —— +X[r , ] —— + ... + X[r/ . (3.5) 

' 3 ' 2 ' 4 ' 3 1 N I jv — 1 



8 



Recall that v Tj is the Dirac measure at r = rj . Clearly, a Dirac measure v Tj , 
means that there is an infinitely thin shell at r = rj with unit ADM mass 
and we will call such a configuration a Dirac shell. The proof of (|3.4p is a 
consequence of a crucial monotonicity property of F as r' Q — > r\ and r[ — ► r-i- 
The next step in our strategy is to show that 

F r M <FrM), 

where v' 2 is the measure obtained by moving the Dirac shell at r = r\ to 
r = r 2 , i.e., 

v'2 = (cm + C2r 2 )4 2 + U2(lr. 
It will be seen that the structure of F allows one to continue this process so 
that the next step is to replace the step function on the interval [r' 2 , 7*3] by 
a Dirac shell with weight 03^ at r = and again show that F applied to 
this measure increases the value. Then we move the Dirac shell with weight 
c\r\ + C2?*2 at r = r2 to r = r% and thus obtain a Dirac shell at r = r% with 
weight c\r\ + C2r2 + 037*3. This measure thus takes the form 

v 'z = ( c i r i + c 2^2 + c 3 r 3 )4- 3 + u 3 dr, 

where 

/ x c 4 r 4 c 5 r 5 c N r N 

Mr) = Xlr> 3 ,r 4 ]Z— 7 + X K ,r s ]—— 7 + - + 



— r' 3 ~ l v ""r^ — r\ ' l ' "~^ TN — r' N _ x 



In this way we obtain the chain of inequalities 

^.(/lu) < FrM) < F rM) <~<Fr*{v' N ), 



where 



Now 



N 



v' N = ^Cjrj 4-, =: m*4„,- (3.6) 

3=1 



F (1/ ) = ~ 2m */ r * ^ 3 7 ^ 

1 + s/l — 2m,/r» 

which follows by using the method in [3], and also from the proof given in 
section 5. In view of (12.80 we thus obtain 



2(1 + 0)771*^1 -2m*/r* , 00 . 

777* < , = , (3.8) 



1 + - 2m*/r 

and solving for 2m* /r* gives 



2m. (1 + 2^-1 
r* < (1 + 20)2 • ^ 



9 



4 An electrostatic analogy 



A classical problem in electrostatics is the question how a unit amount of 
charge should be spread over a bounded set £ G I 3 in order to minimize 
the Coulomb energy 

£{p)-=\( I P{x)p{y)\x - y\ 2 ~ n dxdy. 



2 



E JE 



Following the exposition in [20] the minimum energy is defined to be 
iCap(-E) -1 , where Cap(-E) is the capacity of E, i.e., 



2Cap(£) 



M{S(p): [ p = l}. (4.1) 
Je 



A minimizing p does exist if E is a closed set. It is not a function but 
a measure (an equilibrium measure) concentrated on the surface of E. In 
particular, if E is a ball or a sphere of radius R then the optimum distribution 
for the charge will be 

' (4-2) 



4vr^ 2 
and 

Cap(B fl ) = R. (4.3) 

Of course, this problem can equivalently be formulated as a variational prob- 
lem for Newtonian gravity but since we wish to stress the relation to capacity 
theory which originates from the electrostatic problem we have preferred to 
use that formulation. 

The analogy with our case should be clear in view of (14. 2p . Let us also 
note that capacity can equivalently be defined as the largest charge that can 
be carried by a body (e.g. a ball with radius R) if the voltage drops by at 
most one, cf. [2]. This formulation suggests that we in our situation define 
the capacity of a ball with radius R to be the largest ADM mass that a 
spherically symmetric static body with area radius R can have. Using this 
definition we then get in view of Theorem 2 that the capacity is given by 

({l + 2fl) 2 -1)R 
2(l + 20) 2 ' 

Of course, we could also introduce a similar definition as in (|4.ip by using 
a variational formulation for F instead of £ . The following theorem, taken 
from [20] , is an interesting feature of balls in R n for the capacity in (|4.ip . 
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Theorem 3 (}20$ ) Let E C R n , n > 3, be a bounded set with Lebesgue 
measure \E\ and let Be be the ball in R n with the same measure. Then 

Cap(B E ) < Cap(E). 

This theorem suggests that spherical symmetry might be an important case 
also for the compactness ratio "2M/R" (assuming one has a proper definition 
of such a quantity) of more general static objects. 



5 Proofs 

Proof of Theorem 1. Consider any admissible triplet, so that in particular 
< 4irr 2 p 6 Lj oc , and let / := Airr 2 p. These are the only conditions of 
an admissible triplet needed in this section, the remaining conditions have 
already been invoked to derive the relations in section 2. We will show that 
(|3.1|) implies (|3,2|) . Hence, assume that there is a r„ > with the property 
that (|3.ip holds. By continuity we can choose r* so that 2m(r*)/r* is as 
close as we wish to the critical value and we choose r* so that 

(1 + 20)* -l < 2Mr A< l_ (l+ (1 + 2^ -lx ;= Q 



(1 + 20) 2 r* 2\ (1 + 20) 

In what follows we use the notation m* := m(r*). Fix e > 0. Let h be such 
that h = on [0,6) and h = f on [6, r*], 5 > 0. Obviously, for a sufficiently 
small 5 > the difference < rrif{r) — rn^{r) is arbitrary small and since 
the integration interval [0, r*] is finite it holds by a continuity argument that 
there is largest 5 > such that \F rt {h) — F Tt {f)\ < e/2. 

Now, since the operator F consists of a composition of integrations, there 
is a natural number N, a sub-division {ro,ri, ...,r^}, rj = 5 + j(r^ — 5)/N, 
of the interval [S, r*], and positive constants {c\, C2, cjy} such that the 
function h defined by 

T ( s cm c 2 r 2 c N r N 

K r ) =X[r ,r 1 ] l"X[ri,r 2 ] ^ - + X[r N _ u r N ] , (5-2) 

satisfies |-F r „(^.) — F r ,(/i)| < e/2, and — m^(?"*)| < e. Here X5 is the 
characteristic function, i.e., xs( r ) = 1 if r £ £?, and xs( r ) = if r ^ S 1 . 
The condition that the constants Cj are positive is technical and it is not 
required that / must be positive, only non-negative, but since we only seek 
an approximation our positivity condition is easy to satisfy. For technical 
reasons we also require that N is taken large, i.e., 
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We now define 

hf ^ Cin i c 2 r 2 , . c N r N 

h{r) = x ^]^r + XKrt^^r + - + x K _ v r N] rN _ ^ • (5-4) 

Here rj < r^- < Pj+i. Note that h = h\i r 1 - = tj for all j € N. Moreover note 
that 

hdr = / /idr, 
/o Jo 

so that the quasi-local mass at r = r* given by the energy densities p = 
h/(4:7rr 2 ), and /? = /i/(47rr 2 ), are the same. The function h, will be the 
main object below. As explained in section 3 we will modify h, by varying 
the parameters r'j and moving parts of the matter, and finally obtain the 
inequality 

F r * (/) < F r , (h)+e< F r . (u' N ) + e, (5.5) 
where u' N is given by (|3.6|) . The proof is split into four steps. 

5iep i. 

In the first step we will by a straightforward computation find an expression 
for F rt (h). Since this computation is crucial and quite lengthy we will present 
the main steps. In what follows j and k will always be non-negative integers. 
Let Co = 0, and let k > 1. From (|5,4p we get 

j=k—l i i \ 

c k r k\P ~ r k-V 

m(a) = y. c j r j H ', > where r^-i < cr < r^. (5.6) 

j=0 ' " r fc -7 Vi 

By defining 

j=fc 

M fc := ^c,-rj, fc > 1, 
j=o 

we thus get 

c^r^a — rl 1 ) 

m(cr) = M fc _i H ? , r fc _i < a < r k . (5.7) 

Next we define 

r 00 h(a) da 
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Note that this is the main expression in the operator F, cf. equation (|2,7p . 
From (|5.4p it thus follows that for r^_ 1 < r] < rj, j > 1, 

g[h\ (v) = r — ^ — -— -T- + r+i cj+irj+i da 



v (r 3 -rj>(l-^) M { r 3+1 -r>)a(l- 2 -^) 
f TN cntn da ~ , . 

+...+y - — ==G i (i y )+g i+1 +...+g Jy . 

■V^-i ( rjv - r^-Jo^l ^ J 

Here the twiddle over the first term emphasizes that it depends on r\ whereas 
the remaining ones do not. By inserting the expression (|5.7p for m we get 



CjVj da 



V (r- r> \a(~\ 2M ^ 2c i r >- r j-i) 

Cjrj da 

Icjrjr'j^ - 2M i _i(r i - rj_ t ) - (7(2^ - r,- - r^) 



Note that the denominator in the integrand is positive in view of (|5.8p . Let 

A? := r j — r j-i> we then get 



(5 = -qrj b / 2c j r j r j-i - 2M j-i A i ~ r j( 2c j r j ~ A j) 
' 2c,t,- - Aj ° g V 2v , - 2.U ; ,A ; - r ? (2c j r J - Aj) 

Analogously we get for the rj independent terms 



G 



j j 



j 



/ 2 W ^ 1 -2M ? _ 1 A J -r J (2c J r J -A J ) v 
b \2ct,t1- 1 - 2Af„-_i A,- - r'- , f 2c,n - A,l / ' 1 ' ; 



2c jrj - A, V 2r ; ,-,/•; ; - 2M i _ 1 A j - r^_ 1 (2c i r j - A,) 

Let us now consider the operator F. From the expression (|5.4f) we have 

N 



Fr * {h) = A~ , Cin 6XP ( " 61 (r?) " £ ^) dr? 

+ ~K~2 J 1 C2r2eXP ( ~ ^ 2 ^) ~ 5Z G J') dr? 



+ ... 



/ cjyrjvexp G N (rj))dr]. (5.11) 
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Since the only dependence on r/ in the integrand is in Gj we thus obtain 



cmexp ( - Y!j=2 G i 



./-;,(/;)- - 2 3/ I <>xp( (</))'/>/ 



c 2 r 2 exp ( -EiLa^i) p 



A, 



^ exp (- G 2 (77)^77 



+^ f " exp ( - G N (ri))dr,. (5.12) 
The first two terms in this expression can be written as 

Ciri6 ~ G2 r^itodq + Sp. fe-^W^e-^G, (5 . 13) 



Ai J r > A 



2 Jr! 



As explained in section 3 the idea is to show that F rt (h) is dominated by 
F rt {y2)-, where V2 is the measure 

!/ 2 ( r ) = Ciri ^ i+C2r . 2 ^ 2+X[r , r3] + ... + x ; (5.14) 

'3 '2 'iV ' jV-1 

and then to show that F Tt (v2) < -F rt (z/ 2 ) where 

f 2 (0 = (ciri + c 2 r 2 )^ 2 +X[4 > r 3 ]— — 7 + -+X[r^_ 1I r J v]— — 5 ■ (5-15) 

'3 '2 'N ' N-l 

The measure z^ 2 can thus be thought of as a modified h where c\ and C2 
have been replaced by c[ = and c' 2 = (ciri + c 2 r2)/r2 respectively, and 
where the limit r' x — ► r 2 has been carried out. Note that the quasi-local 
mass generated by v' 2 and h are the same, i.e., m^'(r*) = m/ l (r*). In order 
to show that F Tt (h) < F^ly'^), the terms in the bracket in (|5. 13j) must be 
dominated by 



lim 

r'^n A 



/ e~ d ^ ) dr,+ lim ^ / e-^drj), (5.16) 

1 Jr^ r^->r 2 A2 J r ^ ' 



which in turn must be dominated by 



l im f 2 e'^dr). (5.17) 
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Here G'i denotes the G— function which corresponds to the measure v' 2 - The 
structure of F{h) revealed in (|5.12p then shows that this procedure can be 
continued: we define the measures z/3 and f 3 by 

n , , \ ° I o c 4 r 4 c N r N 
u 3 (r) = (c 1 ri + c 2 r 2 )i^+c a r 3 i^+xy a ,r^- ZF + - + X[ r ' N _ 1 ,r N ]- ~, > 

'4 '3 >N ' N-l 

1 ( \ ( , 1 \ o . c 4 r 4 cn^n 

u 3 (r) = {cm + c 2 r 2 + c 3 r 3 )^3 + X[ r ' 3 ,r 4 ]- -7 + - + X[r' N _ v r N ]——i , 

' 4 — '3 * N ~ 1 jv_i 

and we show that F rt {h>'^) < F rt (v 3 ) < F rt {v' 3 ). In this way we obtain a 
chain of inequalities 

FrM < F r M) < FrM) < - < F rM), 

where v' N is the Dirac measure at r = rjy = r* with ?n^(r*) = m/^r*). Let 
us now compute the sum of the two terms in the bracket in (|5. 13j) . We use 
the following notation 

and 

To = 

A 



e -^V)dr), (5.18) 



■/'., = ^ e-^dq. (5.19) 
1 
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We have from (|5,9p 



f ri / 2c iri r' - n(2ciri - AQ x ss^hsj dj? 
7^ V 2cmr^ - 7/(2cm - Ai) / ''' 

2cmr' - ri(2ciri - A x ) 



cm 

2cm-Ai 



(2cm - Ai) 



2cinr' Q - r?( 2c m - Ai) 



2c i r i — ^1 



2cm7^ - ri(2ciri - A x 



c m 



cm 



Ai 



^(2c iri r' - r' (2c iri - Aj 



2cirir' — ri(2em — A\ 
^2cmr' - ri(2cm - Ai) 



2 =n 



cm - Ai 

2arir' Q - r' (2ari - Aj 



{( 



c 1 r 1 -A 1 



^2drir' Q - n(2cin - A^ 

Furthermore, from (|5.10p we have 

e _G 2 = / 2c 2 r 2 r; - 2MjA 2 - r 2 (2c 2 r 2 - A 2 ) 
V2c 2 r 2 r; - 2MiA 2 - ri(2c 2 r 2 - A 2 ; 

The term Ti can thus be written 



.}. (5.20) 



e 2 r 2 



2c 2 r 2 -A 2 



cm 



^2cirirQ - n(2cm 



Ai 



cm 



Ai 



/ 2cirirg - rp(2cm - Ai) 



c 1 r 1 -A 1 
2c 1 r 1 — Aj 



V2cmr^ - ri(2cm - A x )) 1 } 
/ 2c 2 r 2 r^ - 2MxA 2 - r 2 (2c 2 r 2 - A 2 ) \ 2e2 %?A a 
\2c 2 r 2 r^ - 2MiA 2 - ri(2c 2 r 2 - A 2 ) 



(5.21) 
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A very similar calculation shows that 

C2? , 2 \2c 2 r 2 r' x - 2MiA 2 - r 2 (2c 2 r 2 - A 2 



To 



A 2 c 2 r 2 - A 2 

J fO„ „ A \ . 'J'J 



f / 2c 2 r 2 rl - 2MxA 2 - r[(2c 2 r 2 - A 2 ) x 2^p4 o 
U2c 2 r 2 ^-2MiA 2 -r 2 (2c 2 r 2 -A 2 ); J 

The aim is to obtain the inequality F rt (h) < F Tt ,(i/ 2 )- Since /i and v' 2 are 
identical for r > it follows from (|5.12p that it is sufficient to obtain the 
estimate T\ + T 2 < T^ 2 + T 2 2 , where T^ 2 and T 2 2 are the corresponding 

terms for v' 2 . Clearly T" 2 = since c'i = 0, and T 2 2 is the expression (|5.17p 
which in view of (|5.22p and the fact that M\ = in this case since = is 
given by 

T 2 ' 2 = lim T 2 , (5.23) 



where 



v r I 2c 2 r 2^i - r 2 (2c' 2 r 2 - A 

To — 



A 2 c 2 r 2 - A 2 

f/ 2c 2 r 2 rl-rU2c 2 r 2 -A 2 ) xfe^ , 
X U2c 2 r 2 ri-r 2 (2c 2 r 2 -A 2 )J 

Here c' 2 = (c\r\ + c 2 r 2 )/r 2 . The expressions for T\, T 2 and T 2 will now be 
simplified. Let us introduce the notation 

r[ 

b k = — =-, fe = 1,2, ... 

which implies that 

A fc = r fc (l - 

Let us consider the term T 2 . By dividing both the numerator and the de- 
nominator by 2c 2 r|, the first factor in the expression (|5.22p can be written 



Car2 [2c 2 r 2 r[ - 2MiA 2 - r 2 (2c 2 r 2 - A 2 ) 
"A 



■2 




c 2 r 2 - A 2 








c 2 r 2 


(b 2 


Mi(l-6 2 ) (l 

C2T2 V 


I-62 
2c 2 


)) 




l-b 2 




C2V2 l — b'2 
2C2T2 2C2 








c 2 r 2 








c 2 r 2 (l 


-2 Cl ^-2c 2 ) 


l-b 2 




1 l-6 2 

2 2c 2 




c 2 


- (1 - 62) 



. (5.25) 
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The second factor can be simplified in a similar way 



c 2 -(l-h 2 ) 



( 2c 2 r 2 r[ - 2M 1 A 2 - r[(2c 2 r 2 - A 2 ) 
\2c 2 r 2 r[ - 2M 1 A 2 - r 2 (2c 2 r 2 - A 2 ) 

'2-^(l-&2)-(l-^) 



2c 2 

V2c 2 ca-ra^ 1 " 2 A 2c 2 -(i-6 2 ) 

h 2 -2c^ ^ t 



l-2 Cl ^-2c 2 



(5.26) 



In conclusion To can be written 



c 2 r 2 (l - 2ci£ - 2c 2 ) r , 1 - 2 Cl a - z 2 



c 2 -z 2 Ul-2 Cl ^-2c 2 



(5.27) 



where we have introduced the notation 

z k = 1 - b k , k = 1,2, ... 
Simplifying T\ and T' 2 in a similar way leads to the following expressions 

= end -2ci)f gSr _ ,| - »»S - *» \ «fe (Mg) 

ci-zi lVl-2ci/ /Vl-2ci^-z 2 / 



and 



r / = (cm + c 2 r 2 )(l - 2 Cl % - 2c 2 ) ^ l - Z2 x iSjrg, _ i 
2 c 2 + ci^-z 2 Ul-2 Cl ^-2c 2 ; J - 

(5.29) 

Note that the expression for T' 2 is obtained from (|5.27p by putting c\ = 
and replacing c 2 by c 2 +c\ri/r 2 in accordance with the previous discussion. 

Step 2. 

In this step we show that F rii {h) < F r:t (v 2 ) by showing that F is monotone 
as 7q — * 7*1 and — ► r 2 , i.e., as zi — ► and z 2 — > 0. Let us define 
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and 

B(z,c) = {^)^, (5.31) 

where z £ [0,1/10], and c G (0,Q/2). Recall the definition of Q in (jBTTj) . 
These are the fundamental functions in the expressions for T\ and T 2 , namely 

in -*„(!- fcOAta, ( 1 _ 2 ;,. i/r2 . • own) 

and 

T ^ _ C2r 2 (l-2ci^-2c 2 ) ^/ z 2 c 2 \ 

1 — 2ciri/r 2 V 1 — 2ciri/r 2 ' 1 — 2ciri/r 2 / 

Let us now see that the domain of definition of the functions A and B is 
relevant. Since is the smallest r with 2m* /r* = Q we have in view of 
(|5.ip cin + c 2 r 2 < r 2 Q/2, and since Q < 1, Qcin + c 2 r 2 < r 2 Q/2, which 
implies that 

c 2 r 2 < ^2(1 - 2cm /r 2 ), 

and thus 

02/(1 - 2 Cl n/r 2 ) < Q/2. 

Since ci < Q/2 it follows that the second argument in A and B in (|5.32p and 
(|5.33p is less than Q/2, i.e., c S (0, Q/2). To see that the first argument in 
the functions A and B belong to [0, 1/10] we first check that the condition 
(|5.3p implies that 

8 + (r* - <5)/iV 9 + Q 



5 + 2(r* - cf)/iV - 10 
The inequality above can be written 

(8 + 2Q)(r,-<5) < (l-Q)tf 



10N ~ 10 

which clearly is satisfied if 

lOr* 



N > 



(1-Q)8 

In view of (|5.3p we thus have for j > 1, 



r'j ^ Tj ^ n _ 8 + (n - <J)/AT ^ 9 + Q 



> -^£- > — = - - y — > — TT\~ ' ( 5 - 34 ) 



r i+ i ~ r j+1 ~ r 2 8 + 2(r* - S)/N ~ 10 



19 



so that Zk = 1 — r' k _ 1 /rk < (1 — Q)/10 for all fc. It follows that 

Z2 Zo 1 



l-2cin/r 2 l-Q-10' 

which proves our claim, i.e., z G [0, 1/10]. 

It is clear that these facts hold in general, i.e., not only for the terms T\ 
and T2 but at any step in our chain of inequalities since X^=i c j r j < r kQ/2. 
We can of course also express the term T' 2 in a similar way but it is not 
useful here. By construction the functions A and B are continuous in the 
domain of definition, in particular they are continuous along the lines z = c 
and z = 2c. 

Lemma 1 For any c £ (0, Q/2) the functions A(-,c) and B(-,c) are de- 
creasing in z, z G [0, 1/10]. 

Proof of Lemma 1. Monotonicity of A. Let us introduce the new variables 

13= and k= (5.35) 

H c I -2c y ' 

We thus have that 

We now express A in terms of these variables and by abuse of notation we 
denote this function again by A. Since 

1 — z 2c — z 

1 + — tt~ = 1 + k(3, 



I -2c 1 — 2c 

it follows that 



A(J3, k) = ^^{(1 + - l}- (5-36) 

We now want to show that d@A > since d z (3 is negative. A straightforward 
computation gives after some rearrangements 



(l-/3) 2 /3 2 (l + fe/3) 
+(J3 - 1)(1 + fc/3) log (1 + + (1 - /3) 2 fc/3] . (5.37) 

Let us denote the factor in square brackets by Adding the first and the 
last term in this expression gives 

V(P,k) = /3 2 (l+A:/3)^ + (/3-l)(l + fc/3) log (1 + fc/3)-/? 2 -fc/3(2/3-l). (5.38) 
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Let 

log (1 + *;/?) ( , 
7 = p , (5-39) 

which is well defined also when (3 = since lim^^oT = k. Since 

2c- z 



k/3 



1 -2c 7 



it follows that k/3 < 1/(1 — Q), and since k/3 is positive as long as 2c > 1/10 
a rough estimate gives 

k/3 > -1/10, (5.40) 

by the condition that z < 1/10. We will below distinguish between the two 
cases < f3 < 2, and (3 < 0. In both cases 7 > 0, or more precisely, in the 
former case we have 7 £ [log (1 + 2k) /2, k], and in the latter case 7 £ (0, k]. 
By using the relation 

k/3 = e 7/3 - 1, 

takes the form 

¥(0,7) = -1 + 2/3 + /3 2 (e 7 -l)+ e 7/3 [(/?-l)/?7- 2/3 + 1]. (5.41) 

By expanding the exponential functions using the formula e x = 1 + x/1! + 
£ 2 /2! + ... and collecting the terms corresponding to different powers in 7 
gives 

•Cft 7) = ^ £ [i-^(^-^ ( u^ji - 1)] V. (5.42, 

Note that the lower orders of 7 vanish. We denote the factors in square 
brackets by and these can thus be written as 

^09) = i (l - (j - \W- 2 + (j - 2)/3^ 1 ) . (5.43) 

We now claim that 

n — /Ti 2 

= -, (1 + 2 ^ + 3/3 2 + + U ~ 2)/3 J - 3 ), j > 3. (5.44) 

This statement is easily shown by an induction argument. First, if j = 3 we 
have from (15.430 that 

<K 3 = |(l-2/? + /3 2 ) = i [ (l-/3) 2 , 
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so the claim is true for j = 3. Assume now that for any positive integer 
P > 2, 

(1 - I3 P {P + 1) + /3 P+1 P) = (1 - /3 2 )(1 + 2/3 + 3/? 2 + ... + P^' 1 ). (5.45) 

We then have by (|5.45[) 

(l-/3 p + 1 (P + 2) + /? p+2 (P + l)) 

= (1 - /3 2 )(1 + 2/3 + 3/? 2 + ... + P^- 1 ) 

+(3 P (P + 1) - 2(P + l)/? p+1 + (3 P+2 (P + 1) 
= (1 - /3 2 )(1 + 2/3 + 3/3 2 + ... + PP P ~ X ) + /3 P (P + 1)(1 - /3) 2 
= (1 - /3 2 )(1 + 2/3 + 3(3 2 + ... + (P + 1)/3 P ), 

and the claim (|5.44p follows. In conclusion we have shown 

oo 

d p A = (1 + kp)ir ^[1 + 2/3 + 3/3 2 + ... + (j - 2)/3^V'- (5-46) 
i=3 J ' 

Note that the lower orders of 7 have vanished. Now, 1 + fc/3 > 0, since 
fc/3 > —1/10, and 7 > 0, so in the case (3 > 0, it follows immediately that 
<9g^4 > 0. Let us therefore consider the remaining case /3 < 0. First we note 
that (3 < implies that z > 2c. Now, since 2; < 1/10 this means that (3 is 
only negative if c is small, i.e., c < 1/20. Therefore, since 7 < k we get 

7 < fc = < 1/18. (5.47) 
1 — 2c 

From the inequality (cf. [T]) 

|log(l-x)| < — , 0<x<l/2, 

we have 

| 7 /3| = I log (1 - fc|/3|)| < 3*1/31/2 < 3/20, (5.48) 

where the last inequality followed from (|5.40p . Let us now estimate the sum 
in (|5.46p . For this we use that 

1 7 7 2 1 „ 2 .11 1 

I + 5l + ¥ + - < 4! (1+ ^ + 7 + - ) = I(T^) < 20' 

by (|5.47|) . together with the formula 

9 1 
1 + 2z + 3aT + ... = -77, -1 < x < 1. 



{l-xf 
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We drop the non-negative terms except for the first one and obtain 

^ + i(l + 2/?) 7 + 1(1 + 2/3 + 3/3 2 ) 7 2 + ^(1 + 2/3 + 3/3 2 + 4/3 3 ) 7 3 + ... 



1 

3! 
1 



4(7|/?l) ; 



1 



7 , 7 



6! + 7! + 8! 



1 



3! 20 L 
1 1 

3! ~ 20 L 
1 



2 7 |/3|+4( 7 |/3|) 3 + 6( 7 |/3|) 5 ...] 



1 + 2 7 |/3|+3( 7 |/3|) 2 +4( 7 |/3|) 3 + ...] 



1 



> 



20 2 



3! 20(1- 7 |/3|) 2 " 3! 20-17 2 



> 0. 



(5.49) 



In the second last inequality we used (|5.48p and in the last (|5.40p . Thus 
dpA > also in the case when (3 < 0, and the monotonicity of A(-,c) 
follows. Let us now turn to the monotonicity of B(-,c). 



Monotonicity of B. We express B in the variables k and (3 and, by abuse 
of notation, get 

B(/3, k) = (l + fc/Jpr. 

As in the case of the function A the claimed monotonicity follows if we can 
show 

d p B > 0. 

We have 



dpB = p±M^ [(l + k 0) log (1 + k(5) - kp]. 

Using the variable 7 defined in (|5,39p we can write the factor in square 
brackets as 

[(1 + kf3) log (1 + fc/3) - k/3] = 7 /3e 7/? - (e 7/3 - 1). 

By letting a = 7 /3 we have a function of one variable and it is elementary 
to show the non-negativity of this expression for any a. This completes the 
proof of the lemma. 



□ 



Step 3. 

In this step we show that F(u2) < F{v' 2 ). Hence, we want to show that 
< lim T' 2 - lim ( lim tA - lim T 2 =: f' 2 -f x - f 2 . 



22^0 V Zl^O 



22^0 
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We have from (IQ7D - (|Q9|) 

fx = r 1 



1 - 2ci 



1 - 2ci 



l-2ciri/r 2 - 2ciri/r 2 



1 - 2c', 



2' 



T 2 = r 2 VI - 2cir!/r 2 - Vl - 2c' 2 Ul- 2c 2 , 



and 



Hence 



r 2 ' = r 2 (i 
f' 2 - % - f 2 = 



l-2c' )\ \-2cL 



1 - 2c' 



r 2 (l - a/ 1 - 2ciri/r 2 ) 



1 - 2ci 



1 - 2ciri/r 2 



(1- vT=2^) 



Define 



then 



k := 1 , so that k G (0, 1), 

^2 



Vl - 2cin/r 2 = V1-2cia/1 + 



1 - 2ci 

The factor in square brackets above can be written 



T: = r 2 1- Jl-2ci 



T2 



l-2ci 



(1- VT 3 ^") 



f'2 



(l- VT^TVl + y 



2cik 
~2^ 



-ri 



1 + T 



1C\K 

-2 C1 



:(1 - V 7 !^ 7 ^) 



?'2 



2ci« 
-2ci 



1 + 



2cik 
1 - 2ci 



VI - 2ci 1 + 



2ci« 
1 - 2ci 



-(l-K)(l-x/r^2HT) 

r 2 



-i _| 2cik L 
1 l-2ci 



1 + 



2cik 



1 - 2ci VI - 2ci 



+ K - 1 



r 2 



v /1-2ci(1-k) 



(5.50) 



V 7 ! - 2ci(l - k) - k - VI -2ci(l - k)) 
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Let us introduce the notation 



T(c u k) : = ^l-2 Cl (l 



so that 



T 2 -T l -T 2 = r 2 \ 



k - VI - 2ci(l - «)), 

-r(ci, «). 



(5.51) 



1 - 2c' 2 



1 -2ci 1 - k) 



We want to show that the right hand side is non-negative for any admitted 
choice of the parameters c\, c 2 and k. Since r(0, k) = the statement follows 



since d Cl T > 0. Indeed, we have 



dci 



(1-zc) 



VI - 2ci ^1 - 2ci(l - re) 
which is positive since re G (0, 1). Hence — T\ — T 2 > 0. 



Step 4- 

At this stage it is clear that by repeating the arguments we obtain 

FrM < KM) < FrM) < - < ^r*Kr), 

where v' N is the Dirac measure at r = = with m v > (r*) = mh{r*). An 
appropriate method for computing F Tt (v' N ) is given in [3]. However, we can 
also use the formula (|5.29p with c\ = 0, c 2 = m/V*)/ 7 "* and z 2 = 0, and we 
get with := m/j(r*) 



FrM) = r.(l - gM-r= - '} - — V = £• 

/-^ 2m, J ^ , , '2m 



The inequalities f|2.8j) and (|5.5p then gives 



2(l + fl K^ 

i + V 1 - 2m '*/ r * 



Using that |m* — m*| < e we get 



2(1 + »)m < . v /l - 2m,/r t 

m* < , = h 0(e), (5.54) 

1 + v 1 — 2m* /r* 
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and solving for 2m»/r» gives 

2m* (l + 20) 2 -l 

< ft + gy +0(e) - (5 ' 55) 

Since e > is arbitrary this contradicts our assumption on 2m*/r*, which 
completes the proof of Theorem 1. 

□ 

Proof of Theorem 2. The proof is a direct consequence of the discussion 
leading to (|2,9|) and the formula ([5.52j) . cf. also [3J. Indeed, let 

2M _ (1 + 20) 2 - 1 
~R ~ (1 + 20) 2 ' 

The formula (|5,52p with r* = R and m* = M gives 

^K) = ™ - , = tt^t, (5.56) 



□ 



1+2Q 

and the proof of Theorem 2 is complete. 
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